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Velocity-Dependent Eddington Factor in Relativistic Radiative Flow 
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Abstract 

We propose a variable Eddington factor, depending on the flow velocity v, for the relativistic radiative 
flow, whose velocity becomes of the order of the speed of light. When the gaseous flow is radiatively 
accelerated up to the relativistic regime, the velocity gradient becomes very large in the direction of the 
flow. As a result, the radiative diffusion may become anisotropic in the comoving frame of the gas. Hence, 
in a flow that is accelerated from subrelativistic to relativistic regimes, the Eddington factor should be 
different from 1/3 even in the diffusion limit. As a simple form, the velocity-dependent Eddington factor 
may be written as f{(3) = 1/3+ (2/3)/?, where (3 = v/c. Using the velocity-dependent Eddington factor, we 
can solve the rigorous equations of the relativistic radiative flow accelerated up to the relativistic speed. 
We also propose a generalized form for a variable Eddington factor as a function of the optical depth r as 
well as the flow velocity: /(r,/?) = 1/3 + (2/3)[l + (t + l)/3]/(l + r -f /?) for a spherically symmetric case. 
The velocity-dependent Eddington factor can be used in various relativistic radiatively-driven flows, such 
as black-hole accretion flows, relativistic astrophysical jets and outflows, and relativistic explosions like 
gamma-ray bursts. 

Key words: accretion, accretion disks — astrophysical jets — gamma-ray bursts — radiative transfer 
— relativity 



1. Velocity-Dependent Eddington Factor 

In the fundamentals of radiation transfer in the static 
atmosphere, in order to solve moment equations truncated 
at the finite order, we need some relation to close the 
sequence (Chandrasekhar 1960; Mihalas 1970; Rybicki, 
Lightman 1979; Mihalas, Mihalas 1984, Shu 1991). When 
moment equations are truncated at the second order, we 
usually adopt the Eddington approximation in the inertial 
frame (laboratory frame) as a closure relation, 
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where E and P'-' are the radiation energy density and the 
radiation stress in the inertial frame, respectively. In a rel- 
ativistic radiative flow, where the gas moves at a relativis- 
tic speed, we also assume the Eddington approximation in 
the comoving frame (fluid frame), 



(2) 



where Eq and Pq are the radiation energy density and 
the radiation stress in the comoving frame, respectively. 
Then, this closure relation in the comoving frame is trans- 
formed to the relation in the inertial frame (Lindquist 
1966; Hsieh, Spiegel 1976; Fukue et al. 1985; Kato et al. 
1998). 

These relations are valid in the diffusion limit, where 
the photon mean-free path is sufficiently smaller than the 
scalelength of the system in the optically thick regime. 
In other words, in such a diffusion limit, the radiative 



diffusion is supposed to be isotropic, and the so-called 
Eddington factor f is set to be 

1 

3' 

However, in the optically thin regime, where the photon 
mean-free path is sufficiently larger than the scalelength 
of the system, the radiative diffusion generally becomes 
anisotropic. For example, in the spherically symmetric 
case, the Eddington factor should be unity in the stream- 
ing limit (it should be noted that in the plane-parallel 
case the Eddington factor is always 1/3). Hence, in order 
to bridge the optically thick to thin regimes, a variable 
Eddington factor, which depends on the optical depth r, 
is usually adopted, 

in the spherically symmetric case (Tamazawa et al. 1975). 

Recently, relativistic radiative flows were rigorously 
solved under a simple situation by Fukue (2005b). He 
examined a flow radiatively-driven perpendicular to a lu- 
minous disk under a fully special relativistic treatment, 
taking into account radiation transfer. He adopted the 
current formalism as well as the Eddington approximation 
in the comoving frame (e.g., Kato et al. 1998). When he 
arranged the basic equations of relativistic radiation hy- 
drodynamics for the one-dimensional flow without grav- 
ity or other forces, he found the singularity in the sense 
of "sonic" points, where the flow speed is equal to the 
relativistic sound speed of c/VS. This singularity origi- 
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nates from the closure relation (2) for the radiation fields, 
and the appearance of the singularity suggest that the 
Eddington approximation would be invalid in the rela- 
tivistic flow with strong velocity gradients (cf. Turolla, 
NobiH 1988; DuUemond 1999). 

In analogy to the optically thick-thin transition, the 
possible invalidity of the Eddington approximation in such 
a relativistic flow can be understood as follows. When the 
radiative flow is accelerated up to the relativistic regime 
and the flow velocity becomes of the order of the speed of 
light, the velocity gradient also becomes very large in the 
direction of flow. As a result, even in the comoving frame 
of the gas, the velocity fields as well as density distribu- 
tions are no longer uniform, and the photon mean-free 
path would be longer in the downstream direction than 
in the upstream or other directions. Hence, the radiative 
diffusion becomes anisotropic, and we should consider a 
variable Eddington factor: the Eddington factor may de- 
pend on the flow velocity v {— f3c) in the form such as 



(5) 



In several literatures (e.g., Mihalas, Mihalas 1984; 
Mihalas 1986; Kato et al. 1998), the influence of the flow 
velocity on radiation transfer was discussed in the con- 
text of the description between the comoving and inertial 
frames. That is, at the lowest order of the flow velocity 
to the speed of light, the closure relation in the inertial 
frame must be modified to the order of (v/c)^. In the co- 
moving frame, however, the Eddington approximation (2) 
was adopted, and no one considered the possible violation 
of the Eddington approximation in the comoving frame. 

In this paper, we thus propose a velocity-dependent vari- 
able Eddington factor for the relativistic radiative flow, 
whose velocity becomes of the order of the speed of light. 

In the next section we describe the basic equations 
for the one-dimensional vertical radiative flow, in order 
to show the typical problem which needs the velocity- 
dependent Eddington factor. In section 3, we exam- 
ine several properties required for the velocity-dependent 
Eddington factor. In section 4, we solve a relativistic ra- 
diative flow under the appropriate boundary conditions 
at the flow base and top, and display the effect of the 
velocity-dependent Eddington factor. The final section 
is devoted to the discussion on the physical reason of 
the present problem, the generalization of the variable 
Eddington factor, and concluding remarks. 

2. Basic Equations 

In order to demonstrate the present problem inherent 
in the relativistic radiation transfer flow, we examine a 
simple one-dimensional flow in what follows. 

We consider the plane-parallel case in the vertical direc- 
tion; a luminous flat disk or thin skin above a luminous 
sphere. The radiative energy is transported in the verti- 
cal direction, and the gas itself also moves in the vertical 
direction by the action of radiation pressure. For simplic- 
ity, in the present paper, the radiation field is sufficiently 



intense that both the gravitational field, e.g., of the cen- 
tral object, and the gas pressure are ignored. The internal 
heating is also ignored. As for the order of the flow veloc- 
ity V, we consider the fully relativistic regime, where the 
terms are retained up to the second order of (v/c). 

Under these assumptions, the radiation hydrodynamic 
equations for steady vertical (z) flows are described as 
follows (Kato et al. 1998; Fukue 2005b). 

The continuity equation is 



pcu ~ J (— const.). 



(6) 



where p is the proper gas density, u the vertical four ve- 
locity, J the mass-loss rate per unit area, and c the speed 
of light. The four velocity u is related to the proper three 
velocity v hy u ~ "fv/c, where 7 is the Lorentz factor, 

7 = VT+^'= \l^/\-{vlcf. 

The equation of motion is 

c^J^ = 1^^1±^ \F^(X + 2«2) _ c{E + P)72«] , (7) 
dz c 

where Kabs and Ksca are the absorption and scattering 
opacities (gray), defined in the comoving frame, E the 
radiation energy density, F the radiative flux, and P the 
radiation pressure observed in the inertial frame. The 
first term in the square bracket on the right-hand side of 
equation (7) means the radiatively-driven force, which is 
modified to the order of u^, whereas the second term is the 
radiation drag force, which is also modified, but roughly 
proportional to the velocity. 

In the no-gas pressure approximation and without heat- 
ing, the energy equation is reduced to a radiative equilib- 
rium relation, 



= j - CKabs£^7^ - CKabs-Pu^ + 2Kabs^'7W, 



(8) 



where j is the emissivity defined in the comoving frame. 
In this equation (8), the third and fourth terms on the 
right-hand side appear in the relativistic regime. 

For radiation fields, the zcroth- moment equation be- 
comes 

dF 

- p-i [j - CKabs-E 4- CKsca{E + P)u'^ 



dz 

+CK^X,sFu/ -f - Ksc!,F{l + / C^)-fu] 

The first-moment equation is 
dP p7 



dz 



-K,cai^(l + 2^2) + CKsUE + P)-fu] 



(9) 



(10) 



In order to close moment equations for radiation fields, 
we need some closure relation. If we adopt the usual 
Eddington approximation. 



(11) 



in the comoving frame as the closure relation (Fukue 
2005b), where Pq and Eq are the quantities in the comov- 
ing frame, the transformed closure relation in the inertial 
frame becomes 



No. 



cP 1 + -u' 
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(see Kato et al. 1998 for details). Instead, we here adopt 
a velocity- dependent variable Eddington factor /(/3), 

Po = fif3)Eo (13) 

in the comoving frame, where /3 = v/c, v being the flow 
velocity. If we adopt this form (13) as the closure relation 
in the comoving frame, the transformed closure relation 
in the inertial frame is 

cP + fu^) = cE {f-/^ - u^) + 2Fyu (1 - /) , (14) 

or equivalently, 

cP (1 - fP^) = cE (/ - p^) + 2F/3 (1 - /) . (15) 

The constraint on the function f{(3) will be discussed later. 

Eliminating j with the help of equations (8), and using 
continuity equation (6), equations (7), (9), and (10) are 
rearranged as 

^du 
dz 
dF 

dz 
dP 



discussed in Fukue (2005b), and for a suitable form of the 
variable Eddington factor /(/3). 

3. Properties of /(/3) 

We here briefly discuss the constraints on the velocity- 
dependent variable Eddington factor /(/3) and choose a 
preferable candidate for /(/3). 

At first, this function must be reduced to 1 /3 in the non- 
relativistic limit of /3 = 0. In the extremely relativistic 
limit of /3 = 1 with strong velocity gradients, where the 
photon mean-free path becomes large, this function would 
approach unity, as in the case of the streaming limit for 
an optically thin regime. Hence, the function /(/?) must 
satisfy the boundary conditions. 
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(25) 
(26) 



CJ— = (Kabs + 'tsca)p- [^(1 + 2u^) ~ c{E + P)ju] , (16) 



(Kabs + Ksca)pw HE + P)jU - F{1 + 2u')] , (17) 



and should monotonically increase from to 1. 

Another constraint comes from the critical conditions. 
As can be seen in equations (23) and (24), these equa- 
tions have singularities (critical points). That is, when 
the condition, 



(Kabs + Ksca)p- [ciE + P)ju -F{1 + 2u')\ . (18) /(/3) -0' = O, 



(27) 



The integration of the sum of equations (16) and (18) 
yields the momentum flux conservation along the flow. 



cJu + P = K (= const. 



(19) 



Similarly, after some manipulations, the integration of the 
sum of equations (16) and (17) gives the energy flux con- 
servation along the flow, 

c^Jj + F = L {= const.). (20) 

At this stage, the basic equations are the equation of 
motion (16), the mass flux (6), the momentum flux (19), 
the energy flux (20), and the closure relation (14). 

Next, by introducing the optical depth t by 

dr = -{K^bs + Ksca,)pdz, (21) 
the equation of motion (16) is rewritten as 

[F{l + 2u^)-c{E + P)-fu\. (22) 



is satified at some f3, the denominator of equation (24) 
vanishes, and the equation becomes singular. For ex- 
ample, if the Eddington factor is fixed as / = 1/3, the 
equation becomes singular at /? = ±l/-\/3, or at the point, 
where the flow velocity is equal to the relativistic sound 
speed of c/\/3- As already pointed out in DuUemond 
(1999), this singularity originates from the finite num- 
ber of moments and the closure relation adopted (see also 
Fukue 2005b). In order for the relativistic transfer flow to 
be always regular, the condition. 



/(/3)-/32>0, 



(28) 



du 7 r , 9s 

dr c 



should be satisfied for 0< (3 <1. Furthermore, since / = 1 
at /3 = 1, the equation is marginally critical {f — 0^ ^ 0) 

■Ai 13^1. 

In addition, from the linear analysis around a critical 
point of equation (23), the velocity gradient near to the 
singularity is found to be 



Furthermore, eliminating E with the help of equation 
(14), this equation (22) can be finally rewritten as 



du 
d7 



du J F(f"f^ + u'^)-cP(l + f)ju 
cJ— = 



^^7^ 1-/ df 



c (1 -/)2u- (l+u2) 



dr c 
or equivalently, 
2 J 2dP 



F{,f + f3^)-cP{l + f)P 



(23) 



(24) 



df 
Uu 



(29) 



where /3 = w/c. When the Eddington factor / is 1/3, equa- 
tions (23) and (24) are reduced to those given in Fukue 
(2005b). 

We then solve equations (23) [or (24)], (19), and (20) for 
appropriate boundary conditions at the moving surface, as 



where the subscript c denotes the values at the critical 
point. If the Eddington factor is fixed as / = 1/3, the ve- 
locity gradient becomes dw/drjc = 3/4 (d/S/drjc = 
Or, the transonic solution for the outflow [u > 0) must 
be decelerated at the critical point, as was discussed in 
Fukue (2005b), where the subsonic solutions were ob- 
tained. Hence, in order for the flow to be accelerated, 
even if the solution is marginally critical (/ — = 0) at 
(3 = 1, the velocity gradient du/drjc should be negative 
(or zero) there for the outflow (u > 0) solutions. Since 
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the numerator on the right-hand side of equation (29) is 
always positive, the velocity gradient is negative if the 
denominator of equation (29) is negative: 



d 



[{i-mi+u')] 



du 

or equivalently, 



d_ 

du 



[(l-/)7^ 



<0, 



d_ 

dp 



1-/ 
1-/32 



<0. 



(30) 



(31) 



Now, bearing these constraints in mind, we exam- 
ine several possible candidates of velocity-dependent 
Eddington factors. 

One simple case is a form of 

f{l3) = l + lf3\ (32) 

where n is constant. In this case, the Eddington factor is a 
monotonically increasing function of /3, which satifies the 
boundary conditions at /3 = and 1. However, since f — 0^ 
becomes negative at some (3 for large values of ri (n > 3), 
the values of n should be restricted as n < 3. Moreover, 
for this form, the equation is marginally critical at /3 = 1, 
and the velocity gradient at the critical point is negative, 
zero, and positive for n = 1, 2, 3, respectively. Hence, the 
case of n = 3 is rejected, and the cases of n = 1 and 2 are 
retained; the simplest case is 71 = 1. 
Alternative is a form of 

m - (33) 

where n is constant. In this case, / is also monotonically 
increasing, but f — 0^ becomes negative for large n (n> 2), 
and the values of n may be around unity. Moreover, for 
this form, the equation is marginally critical at /? = 1, and 
the velocity gradient at the critical point is positive, and 
this form may be inadequate. 

In analogy to a variable Eddington factor for the opti- 
cally thick-thin transition, we can consider a form of 



l + l/'a_ /3+ ^1-/32 



+ 3/u /3 + 3yr3^ 



(34) 



In this case, / is monotonically increasing, but f ~ 0^ 
becomes negative at large u. Similar forms such as / = 
(l-hl/u")/(l-f3/u") or /=[l-hl/(7-l)"]/[l-f 3/(7-1)"] 
would be also inadequate, although we have not checked 
all possible cases one by one. 

Consequently, we adopt the simplest case of 

f{P)^\ + \p (35) 
in this paper. 

4. Examples of Solutions 

In this section we show several examples of solutions 
of equations (24), (19), and (20) with velocity-dependent 
Eddington factors under appropriate boundary conditions 
at the moving surface. 



As already pointed out in Fukue (2005b), the usual 
boundary conditions for the static atmosphere cannot be 
used for the present radiative flow, which moves with ve- 
locity at the order of the speed of light. 

When there is no motion in the gas ("static photo- 
sphere"), the radiation field just above the surface un- 
der the plane-parallel approximation is easily obtained. 
Namely, just above the disk with surface intensity /g, the 
radiation energy density i?s, the radiative flux i^g, and the 
radiation pressure Pg are (2/c)7r/s, Tr/g, and (2/3c)7r/s, 
respectively, where the subscript s denotes the values at 
the disk surface. However, the radiation field just above 
the surface changes when the gas itself does move upward 
("moving photosphere"), since the direciton and inten- 
sity of radiation change due to relativistic aberration and 
Doppler effect (cf. Kato et al. 1998; Fukue 2000). 

Let us suppose a situation that a fiat infinite photo- 
sphere with surface intensity Is in the comoving frame is 
not static, but moving upward with a speed Vs (= C/Sg, 
and the corresponding Lorentz factor is 7g), where the 
subscript s denotes the values at the surface. Then, just 
above the surface, the radiation energy density Eg, the 
radiative flux Pg, and the radiation pressure Pg measured 
in the inertial frame become, respectively. 



cPg = 27r/g 



Pg = 2tiL 



cPg = 27r/g 



37s^ 


+ 37gMg 


+ Ws 




3 




37g^ 


+ 87sWg 


+ 3ug2 
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7s^ + 


37s 4 


3u2 


3 



(36) 
(37) 
(38) 



where Ug (=7gWg/c) is the flow four velocity at the surface 
(Fukue 2005b). 

Thus, we impose the following boundary conditions: At 
the flow base (deep "inside" the atmosphere) with an ar- 
bitrary optical depth tq, the flow velocity u is zero, the 
radiative flux is Fq (which is a measure of the strength of 
radiation field), and the radiation pressure is Pq (which 
connects with the radiation pressure gradient and relates 
to the internal structure), where the subscript denotes 
the values at the fiow base. At the flow top ("surface" 
of the atmosphere) where the optical depth is zero, the 
radiation field should satisfy the values above a moving 
photosphere given by equations (36)"(38). 

Applying boundary conditions (36)~(38) to equations 
(19) and (20), we have two relations on the boundary 
values and mass-loss rate: 



Jc^Us + cPs^cPo, 



Fq. 



(39) 
(40) 



Physically speaking, in the radiative fiow starting from 
the fiow base with an arbitrary optical depth tq , for initial 
values of Fq and Pq at the fiow base, the final values of 
the radiation fields Eg, Eg, Pg, and the flow velocity Wg at 
the flow top can be obtained by solving basic equations. 
Furthermore, the mass-loss rate J is determined as an 
eigenvalue so as to satisfy the bondary condition at the 
flow top (cf. Fukue 2005a in the subrelativistic regime). 
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Fig. 1. Flow three veloeity v (solid eurves), radiative flux F 
(dashed curves), and radiation pressure P (doted curves), as 
a function of the optical depth t for several values of Vs at the 
flow top in a few cases of tq: (a) ro = 1 and (b) ro = 10. The 
values of are 0.1, 0.3, 0.5, 0.7, 0.9, and 0.99 from bottom 
to top of V and 0.1, 0.3, and 0.5 from bottom to top of F and 
P. The velocity v is in units of c, F and cP in tt/s- 

In the present, fully relativistic case, however, the final 
values of the radiation fields at the flow top depend on 
the flow velocity there, and the flnal values at the flow 
top cannot be analytically expressed by the initial values 
at the flow base. Hence, in this paper we determine the 
mass-loss rate as follows. 

In the radiative flow with optical depth tq , we first give 
the final flow velocity (and 7s), instead of the initial 
value of Pq. Then, the final values of radiation fields Eg, 
Fs, and Pg can be fixed by equations (36)-(38). Next, 
we give a trial value for the mass-loss rate J, and the 
initial values of Pq and Fq can be fixed by equations (39) 
and (40). Since all the parameters are temporarily fixed, 
we solve equation (24) from r = tq to t = 0. Generally, 
however, the obtained final velocity at t = is different 
from a given Vs- Thus, we vary the value of J and follow 
iterative processes, so that the calculated final velocity 



Fig. 2. Final velocity Vs at the flow top (solid curves), and 
the mass- loss rate J (chain-dotted ones), as a function of Pq 
for several values of tq at the flow base: tq = 1 and 10. The 
quantities are normalized in units of c and Tr/g. That is, the 
unit of F and cP is tt/b and the unit of J is nls/c^. 

coincides with a given final velocity Vg. 

Finally, we adopt the velocity-dependent Eddington fac- 
tor of the form (35). 

Examples of the relativistic radiative fiows under the 
present velocity-dependent Eddington factor arc shown in 
figures 1 and 2. 

In figure 1 we show the flow three velocity v (solid 
curves), the radiative flux F (dashed curves), and the ra- 
diation pressure P (doted curves), as a function of the 
optical depth r for several values of Vs at the flow top and 
for To = 1 and 10. The quantities are normalized in units 
of c and nig. 

In contrast to Fukue (2005b), where the acceleration 
solution beyond c/ -\/3 was not obtained by the presence of 
the singularity, originating from the usual fixed Eddington 
factor, the present fiow velocity can easily exceed c/V3 to 
be accelerated toward the highly relativistic regime, since 
the singularity is evaded by the adoption of the velocity- 
dependent Eddington factor (35). 

Other properties of the fiow are similar to the previous 
case. For example, when the initial radiative fiux Fq at 
the flow base is large, the fiow is effectively accelerated, 
and the final speed of the flow becomes large. The mass- 
loss rate J is not given arbitrarily, but is determined as 
an eigenvalue (Fukue 2005a, b). 

In flgure 2 we show the final velocity Vg at the fiow 
top (solid curves) and the mass-loss rate J (chain-dotted 
ones), as a function of Pq for several values of tq at the 
fiow base. The quantities are normalized in units of c and 
tt/s. For example, the unit of J is ttIs/c^. 

As can be seen in figure 1, as the radiative fiux increases, 
the final flow velocity at the fiow top increases. Moreover, 
as can be seen in figure 2, and similar to the previous case 
(Fukue 2005a, b), in order for the flow to exist, the radi- 
ation pressure Pq at the flow base is restricted in some 
range. In the subrelativistic case without gravity and 
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Fig. 3. Flow three veloeity v in several forms of veloc- 
ity-dependent Eddington factors. The parameters arc tq = 1 
and Vb = 0.9 c. 

heating (Fukue 2005a), the initial pressure Pq is proved 
to be restricted in the range of 2/3 < cPp/Fg < 2/3 -I- tq. 
In the present case, the initial pressure is also restricted in 
some range, but is modified due to the relativistic effect 
and the velocity-dependent Eddington factor. In addi- 
tion, the loaded mass increases as the initial optical depth 
increases, even if the final flow speeds are the same. 

For comparison, we show in figure 3 the results for sev- 
eral other forms of velocity-dependent Eddington factors. 
In figure 3, the parameters are fixed as tq = 1 and Vs = 0.9 c. 
A solid curve is for the case of /(/3) = 1/3 -t- 2/3/3, whereas 
a dashed one is for the case of /(/3) = 1/3-1-2/3^/3, and 
a dotted one is for the case of /(/3) = 1/(3 — 2/3). The 
boundary conditions at the flow top are the same in each 
case, but the radiation fields {Fq and Pq) at the fiow base 
and the loaded mass J are different in each case. In par- 
ticular, the loaded mass in the case of /(/3) ~ 1/3 + 2/3/3 
is larger than those in other cases examined, and this is 
another reason to prefer the present form. 

It should be finally noted that, thanks to the velocity- 
dependent Eddington factor, the relativistic radiative 
flow can be accelerated to reach an extremely relativis- 
tic regime, and the barrier of the magic speed above the 
static atmosphere (Icke 1989) would be now completely 
cleared away. 

5. Discussion and Concluding Remarks 

Moment equations for relativistic radiation transfer 
have been derived in several literature (Lindquist 1966; 
Anderson, Spiegel 1972; Hsieh, Spiegel 1976; Thorne 1981; 
Udey, Israel 1982; Schweizer 1982). A complete set of 
moment equations for relativistic flow is given by the pro- 
jected, symmetric, trace-free formalism by Thorne (1981). 
Since a moment expansion gives an infinite set of equa- 
tions, one must truncate the expansion at the finite order 
by adopting a suitable closure assumption, in order to 
make the transfer problem tractable. When one truncate 
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at the second order, for example, the Eddington approxi- 
mation is usually adopted as a closure relation. 

Such radiation moment formalism is quite convenient, 
and it is a powerful tool for tackling problems of rela- 
tivistic radiation hydrodynamics. However, its validity is 
never known unless fully angle-dependent radiation trans- 
fer equation is solved. 

Actually, the pathological behavior in relativistic radi- 
ation moment equations has been pointed out (TuroUa, 
Nobili 1988; DuUemond 1999). Namely, moment equa- 
tions for relativistic radiation transfer can have singular 
(critical) points. For example, in the one-dimensional rela- 
tivistic radiation flow using the Eddington approximation 
(2), where the moment equations are truncated at the sec- 
ond order, the singularity appears when the flow velocity 
becomes /3 = u/c= l/V^ (e.g., Fukue 2005b). As a result, 
solutions behave pathologically in regions of strong veloc- 
ity gradients. The appearance of singularities is explained 
because we approximate the full transfer equations with 
a flnite number of moments (DuUemond 1999). 

In order to avoid the singularity, in this paper, we 
thus propose and examine a velocity-dependent variable 
Eddington factor in radiation transfer in relativistic ra- 
diative flows, which are accelerated up to a relativistic 
speed with strong velocity gradients. 

As already stated in section 1, the physical reason is 
clear. In an optically thick, low velocity regime, the mean 
free path of photons are the same in all directions, and 
the radiative diffusion is isotropic. In a relativistically ac- 
celerated flow, however, the velocity gradient of the flow 
becomes very large, and the density distribution also be- 
comes non-uniform even in the comoving frame of the gas. 
Hence, the mean free path becomes longer in the down- 
stream direction than in the upstream and other direc- 
tions, even in an optically thick regime, and the radiative 
diffusion becomes anisotropic. Thus, in the case of sub- 
relativistic to relativistic regimes, as in the case of opti- 
cally thick to thin regimes, we should consider a variable 
Eddington factor which depends on the flow velocity v. 
The most preferable one, we propose, is 

/(/?) = ^ + (41) 

where (3 = v/c. 

By adopting the velocity-dependent Eddington factor, 
the singularity found under the traditional formalism is 
removed, and the relativistic radiative flow can be solved 
from a low velocity regime to an extremely relativistic 
regime. It should be noted, however, that the present 
approach is relevant only for moment expansion up to the 
second order. If we attempt to solve the moment equation 
up to higher order, we may search alternative ways. 

In this paper, we considered only the one-dimensinal 
case. The vector form of the present variable Eddington 
factor may be written as 

when the velocity components are (0,0,/3). 



No. 



Velocity-Dependent Eddington Factor 



We also propose a general form of a variable Eddington 
factor f{T,P), which would be useful for an arbitrary op- 
tical depth r, and for an arbitrary velocity v {= (3c). The 
conditions for such a generalized Eddington factor are: (i) 
in the low velocity limit {(5 ^ Q), it has the similar form 
as equation (4), (ii) in the optically thick limit (r ^ 1), it 
becomes the presently proposed form, (iii) in the optically 
thick, low velocity limit, it is 1/3, and (iv) in the optically 
thin limit (r ~ 0) and/or in the relativistic limit (/3 ~ 1), 
it approaches unity. For a spherically symmetric case, in 
order to satisfy these conditions, we apply the relativistic 
summation rule for the flow speed v and the photon diffu- 
sion speed c/(r + 1), and obtain a generalized Eddington 
factor as 



1 2r- 



1 



T 



T 



21 + {t + \)(3 



3 1 



(43) 



It should be noted that, in the low velocity limit (/3 ^ 
0), this hybrid form (43) becomes /(t,0) = (3-hr)/(3-f 
3t), which is somewhat different with equation (4), but 
essentially the same. In the optically thick limit (t ^ 1), 
this equation (43) is reduced to equation (35). For a plane- 
parallel case, on the other hand, equation (41) can be used 
for an arbitrary optical depth. 

In the present paper, we examine the plane-parallel case 
with the form (41). The spherically symmetric case using 
the hybrid form (43) will be discussed in a seperate paper. 

The velocity- dependent variable Eddington factor for 
relativistic radiative flows proposed in the present paper 
is fundamentally important in various aspects of relativis- 
tic astrophysics with radiation transfer. For example, the 
present form may be applied to the cases of black-hole ac- 
cretion flows with supercritical accretion rates, relativistic 
jets and winds driven by luminous central objects, rela- 
tivistic explosions including gamma-ray bursts, neutrino 
transfers in supernova explosions, and various events oc- 
cured in the proto universe. 

Of course, the present proposition is only the first 
one, and not the final one. There exist many points 
to be improved, including more suitable forms for 
velocity-dependent Eddington factors. For example, the 
Eddington factor may depend not on the velocity, but on 
the velocity gradient. In this field of relativistic radiation 
tranfer, there still remain many problems to be solved. 
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